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Abstract 

BES data on J/tp — > r y(K + K~tt + tt~) are presented. The K*K* contribution 
peaks strongly near threshold. It is fitted with a broad _+ resonance with 
mass M = 1800 ± 100 MeV, width r = 500 ± 200 MeV. A broad 2++ reso- 
nance peaking at 2020 MeV is also required with width ~ 500 MeV. There is 
further evidence for a 2 h component peaking at 2.55 GeV. The non-K*K* 
contribution is close to phase space; it peaks at 2.6 GeV and is very different 
from K*K*. 

PACS numbers: 14.40.Cs, 12.39.Mk, 13.25.Jx, 13.40.Hq 
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In 1990, MARK III presented their spin-parity analysis of J/ip — > ^K*K* at the Rheinfels 
Workshop on the Hadron Mass Spectrum jjfl. They found a dominant h component, 
accounting for 55% of the data. In addition, smaller 2 + and 2~ channels of about equal 
strength were observed. Interferences were not included. Here we present BES data and 
carry out a full partial wave analysis. 

Lattice QCD and other theoretical models predict _+ and 2 ++ glueballs with masses 2.0 
to 2.4 GeV 0]. Recently, several broad structures have been identified in this mass range. 
In an analysis of Mark III data on radiative decays of J/ip to 777r + 7r~, pp, uu, K*K* and (fxf), 
a very broad 0~ component has been found || with a mass of 1750-2190 MeV and a width 
of order 1 GeV. Its decays are flavour blind to first approximation, making it a candidate for 
the 0~ glueball. This 0~ component appears strongly in BES data on J/ip — > 7(7r + 7r~7r + 7r~) 
0| . There, we also find a broad 2 + contribution fitted as a resonance at 1940 ±60 MeV with 
T = 350 ± 100 MeV. The WA91 and WA102 groups have also found a broad 2 + resonance 
at 1920 MeV in their 4n mass spectrum ||. A recent partial wave analysis of Crystal Barrel 
data on pp — > 777771-° |j] has found a broad /2(1980) — > T]i] resonance with mass 1980 MeV and 
width 500 MeV. There are further data on pp — > t]7i°it J7|||, where a broad 0~ component 
is observed plus a peak in rja at 2140 MeV with width ~ 150 MeV. In Refs. || and [f|], an 
/o(2100) was observed in 4ir states. A peak at this mass was observed earlier in the E760 
experiment f9|. The / (2100) decay to rji] has been confirmed in two sets of Crystal Barrel 
data |I|[n] . 



The analysis in this paper uses 7.8 x 10 6 J/ip triggers collected by the Beijing Spec- 
trometer (BES). This detector has been described in detail in Ref ||12|| . Here we describe 



briefly those detector elements crucial to this measurement. Tracking is provided by a 10 
superlayer main drift chamber (MDC). Each superlayer contains four layers of sense wires 
measuring both the position and the ionization energy loss (dE/dx) of charged particles. 



The momentum resolution is ap/P = 1.7%vl + P 2 , where P is the momentum of charged 
tracks in GeV/c. The resolution of the dE/dx measurement is about 9%. This provides 
good it/ K separation and proton identification in the low momentum region. An array of 48 



scintillation counters surrounding the MDC measures the time-of-flight (TOF) of charged 
tracks with a resolution of 330ps for hadrons. Outside the TOF system is an electromagnetic 
calorimeter composed of streamer tubes and lead sheets with a z positional resolution of 
4 cm. The energy resolution scales as a E /E = 22%/y/E, where E is the energy in GeV. 
Outside the shower counter is a solenoidal magnet producing a 0.4 Tesla magnetic field. 

Each candidate event is required to have exactly four charged tracks. Every track must 
have a good helix fit in the polar angle range —0.8 < cos# < 0.8 and a transverse momentum 
> 60 MeV/c. A vertex is required within an interaction region ±20 cm longitudinally and 
2 cm radially. At least one reconstructed 7 is required in the barrel shower counter. A 
minimum energy cut of 80 MeV is imposed on the photons. Showers associated with charged 
tracks are also removed. 

A positive identification of at least one K ± and one 7r ± is required using time of flight 
and/or dE/dx. If two tracks are ambiguous, both alternative identifications are tried. Events 
are fitted kinematically to the 4C hypothesis J/ip — > ^(K + K^tt + tt^), requiring a confidence 
level > 5%. If there is more than one photon, the fit is repeated using all permutations. 
Events with two or more photons are also fitted to J/ip — > 77.fr + K~it + it~ . Those giv- 
ing a better fit than to 'y(K + K~n + iT~) are rejected, as are events fitting the final state 
K + K^7i + ti^ . Next, we require | U m i SS \ = \ E miss — P m i SS |< 0.12 GeV/c 2 , so as to reject 
the events with multi-photons or more or less than two charged kaons; here, E miss and 
P m iss are, respectively, the missing energy and missing momentum of all charged particles; 
they are calculated by assuming the charged particles are K + K~7t + tt~ . The momentum of 
the K + K~ir + n~ system transverse to the photon Pf = 4 | P miss | 2 sin 2 (6> m7 /2) < 0.006 
(GeV/c) 2 is required in order to remove the background J/ip — > ir°K + K~ir + iT~; here 9 my is 
the angle between the missing momentum and the photon direction. Finally, to remove a 
small background of K®, we perform a cut on the ir + n~ invariant mass, | M n + n - —M K o |> 25 
MeV; to remove the small background due to 0(1020), we use a cut on the K + K~ invariant 
mass, I Mk+k- — \> 20 MeV; background J/ip — > ujK + K~ events are eliminated by 
the cut I Mtt+tt-ttO — M w |> 25 MeV, in fitting the ir°K + K~7r + n~ hypothesis with only one 
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photon detected and the 7r° associated to the missing momentum. The number of surviving 
events is 1516. 

The effects of the various selection cuts on the data is simulated with a full Monte 
Carlo of the BES detector; 500,000 Monte Carlo events are successfully fitted to J/ip — > 
, y(K + K~7T + 7T~); all background reactions are similarly fitted to this channel. The estimated 
background is 16%, purely from J/ip — > 7t (K + K~7T + 7T~). We have included this background 
in the amplitude analysis, but it lies very close to non-i^*^* events and has negligible effect 
on the analysis of K*K* events. 

Fig. [l] shows the KK, tctc, Kit, KKtc, K*K and K*tc invariant mass distributions of 
selected events after acceptance cuts. These cuts are accomodated in the maximum likeli- 
hood fit, using the full Monte Carlo simulation of the detector. A strong K* signal is clearly 
visible. There is just a hint of possible structure in Fig. 0(a) around 1500 MeV, but includ- 
ing /o(1500) or /2(1525) into the fit fails to reveal any significant contribution. It is too 
narrow for /2(1525). Figs. 0(e) and (f) show no significant structure attributable to heavier 
Ki or i^*'s states. The K + K~ir + 7t~ invariant mass distribution of J/ip — > r yK + K~ix + it~ is 
shown in Fig. 0(a). The K*K* mass spectrum is obtained with a double cut on the two Kir 
invariant masses, | Mx n — M^* \ < 75 MeV and is shown in Fig. 0(b). Fig. 0(a) shows the 
K + K~tt + 7t~ mass spectrum of J/ip — > K~7t + tt~ when one Kir or both Kit masses fall 
outside the K* region. 

We first attempted to fit the non-K*K* data of Fig. 0(a) as kk and K*°k + c.c; here 



k means the Kir S-wave. The k is parametrized to fit LASS data [13j on ^(1430). Each 
k is made from K + tt~ or K~tt + , and if the k really dominates this process, the Kir mass 
distribution should follow it. The dashed curve in Fig. 0(b) shows that this fails. 

A simple solution which works well is to fit instead with K + K^7t + tt~ phase space. This 
is shown by the full curve of Fig. ^(b). A detailed study of non-i^*^* events reveals only a 
weak p(770) signal, visible in Fig. 0(b). However, when we perform a cut on the 7r + 7r~ mass: 
\M nn — M p \ < 100 MeV, we have been unable to correlate this signal with any particular 
K + K~, K ± tt + tt^ or K + K~7t + tt~ resonance. 



We now turn to K*K* events. One sees immediately a large difference between K* K* 
and non-K*K* events from Fig. |2|(b) and 0(a). The K*K* events peak strongly close to 
threshold. We shall show that the data require a strong 0~ peak at threshold. 

The amplitudes in the PWA analysis are constructed from Lorentz-invariant combina- 
tions of the 4- vectors and the photon polarization for J/ip initial states with helicity ±1. 
Cross sections are summed over photon polarizations. The relative magnitudes and phases 
of the amplitudes are determined by a maximum likelihood fit. We use i to denote the or- 
bital angular momentum between the photon and KKtttt states in the production reaction. 
Because this is an electromagnetic transition, the same phase is used for amplitudes with 
different I but otherwise the same final state. In fitting K*K*, spin-parity assignments up 
to J = 4 have been tried. 

We have examined slices of the K*K* mass spectrum ~ 200 MeV wide. Each slice has 
been fitted with a constant contribution with quantum numbers ++ , 0~ + , 1 ++ , 2 ++ , 2~ + 
and 4 ++ . These different quantum numbers give angular distributions which are distinctively 
different. They depend on two angles (a) 9 K + or 9 n + of K + or 7r + with respect to the K + n~ 
or K~7i + pair in their rest frame, and (b) the azimuthal angle x between the planes of K + tt~ 
and K~7i + in the rest frame of the resonance. Results of the slice fit are shown in Fig. |4[ 
The h contribution is largest and peaks towards the low mass end. There is evidence 
for a broad 2 ++ contribution peaking at ~ 2050 MeV, and this will be confirmed below by 
detailed fits. Contributions from ++ and 4 ++ are small or absent. Though the results of 
the slice fit show that there may be some 1 ++ contributions in the 2.0 ~ 2.4 GeV mass 
region and some 2 _+ contributions at low mass, there is some degree of cross-talk with 0~ 
and 2 + signals; when we put 1 ++ or 2 h into the full fit in these regions, their contributions 
are found to be negligible and are absorbed into 0~ or 2 + . At high masses, there is some 
evidence for a 2~ + contribution. 

The precise mass dependence of each contribution is hard to establish because signals are 
broad and are affected by systematic uncertainties as follows, (a) There is an unknown form 
factor for the vertex J/ip — > 7X. For consistency with other work, we adopt a very weak 
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form factor exp(— ct g q 2 ) with a g = 0.13 GeV -2 , fitted to many channels of J/ip radiative 
decay |TJj]; q is the photon momentum in the J/tp rest frame, (b) Likewise we use a form 
factor exp(—adp 2 ) with ad = 2.0 GeV -2 , where p is the momentum of the final K* in the 
rest frame of X; (c) The resonance X may couple to further channels at higher mass; this 
would pull the cross section down when these channels open. The effect of these form factors 
is to allow shifts of ±50 MeV in the mass fitted to the broad component; however, fitted 
intensities hardly change in shape at all. 

We begin by fitting simple Breit-Wigner resonances of constant width to channels 0~ + , 
2 ++ and 2~ + , including the form factors discussed above. Appropriate Blatt-Weisskopf 
centrifugal barrier factors are included with a radius consistent with the a parameters. The 
best fit is shown on Fig. |^. Crosses are data and histograms the fit. 

The optimum mass for the broad 0~ contribution is 1800 ± 100 MeV. The width opti- 
mises at 500 ± 200 MeV. The present fit with a simple Breit-Wigner form is likely to be an 
over-simplification for such a broad state. In Ref. [14]], all radiative decays are fitted simul- 
taneously. This produces a slightly broader peak, shown by the dashed curve in Fig. |^(b), 
but with little effect on other components. 

The 2 + contribution of Fig. |5|(c) is definitely required. Without it, log likelihood gets 
worse by 58; this is roughly an 8<r effect, from tables of log likelihood for 6 fitted parameters. 
The 2~ contribution is less secure. It improves log likelihood by 18.9 (~ 4.7a). Both the 2~ 
and 0~ contributions are produced in a P-wave and consequently vary as g 3 , where q is the 
momentum of the photon; this dependence leads to a slight underfit at the highest masses. 
If this defect is real, it is most likely fitted by a small El cross section for production of high 
mass + or 2 + components, which cannot be identified definitively. 

The broad 2 + component near 2 GeV is interesting. It is very close to that required to 
fit BES data on J/ip -> 7(7r + 7r~7r + 7r~). It peaks at 2020 MeV. Because the K*K* phase 
space opens rapidly, the pole position tends to come out lower, around 1900 MeV with a 
width of 500 MeV. Changing the form factors moves the pole position up or down slightly, 
but leaves the fitted intensity almost unchanged. 



Figs. H|(a) and (b) shows fitted angular distributions, summed over all K*K* masses. 
Figs. |6](c) and (d) illustrate angular distributions from individual quantum numbers. The 
necessity for a large 0~ contribution is obvious. The full angular dependence involves cor- 
relations between 9 and x (defined earlier), containing much more information than the 
projections of Fig. |^. 

No significant ++ , 1 ++ or 4 ++ resonance is found in the full fit (changes in log liklihood 
< 2.5). That result is itself interesting. There is no evidence for the /q(2100), which appears 
strongly in J/ip — > 7(471") and also in two sets of Crystal Barrel data in the rji] channel and 
also in E760 data. Its absence from the present data suggests it is not an ss state. 

Table | summarises branching ratios evaluated from the fit; the first error is statistical 
and the second covers systematic errors in the overall normalisation of the number of J/ip 
interactions. The r] c component will be the subject of a separate publication, and has been 
subtracted. We find 1516 KKtttt events (of which 320 are K*K*) within the cuts we have 
applied. For each K*, 17.5% of events lie outside cuts. We correct branching ratios to allow 
for this loss. We also correct for all charges in K*K* channels: we have multipied values 
obtained here by a factor (1.5 2 ) for Clebsch-Gordan coefficients in K* decay and a further 
factor 2 for production of K* + K*~ in the radiative decay process. 

We have also searched for a 2~ contribution. A component with these quantum numbers 



was found in J/ijj — > ^(rjir+ir ) at 1840 MeV ||L5 |. A possible interpretation is as an / = 



hybrid; an ssg hybrid is then to be expected around 2.1 GeV with roughly half the branching 
ratio, summed over all kaonic channels. We find that such a 2 _ component with width 250- 
400 MeV produces only an insignificant improvement in the fit: an improvement in log 
likelihood of 4.5 for 4 extra parameters. The branching ratio for production and decay of 
the fitted 2~ component, corrected for all charge states, is then 1.8 x 10 -4 . This is to be 
compared with the value (9.3 ± 3.33) x 10~ 4 observed in all rjinT channels. 

We summarize as follows. We fit with a broad 0~ + resonance with M = 1800 ± 100 
MeV and T = 500 ± 200 MeV, decaying to K*K*. It is to be identified with the broad 0" 
component found also in i]tt + it~ , pp, uouo and 00. The data also definitely required a broad 

8 



2 ++ resonance / 2 (1950) decaying to K*K*. The fit to the high mass end of the K + K~-k + -k~ 
spectrum is best achieved with an additional 2~ + signal. 

The BES group thanks the staff of IHEP for technical support in running the experiment. 
This work is supported in part by China Postdoctoral Science Foundation and National 
Natural Science Foundation of China under contract Nos. 19991480, 19825116 and 19605007; 
and by the Chinese Academy of Sciences under contract No. KJ 95T-03(IHEP). We also 
acknowledge financial support from the Royal Society for collaboration between Chinese and 
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FIG. 1. The KK, tttt, Ktt (two entries/event), KKtt (two entries/event), K*K and K*tt masses 
of J/ip -> -/(K+K-tt+tt'). 
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FIG. 2. (a) The KKtttt mass of J/ip -> 7 K + ii:~7r + 7r'; (b) The K*K* mass of J/ip 7^*^* 




m i?KK + K (GeV/c 2 ) (GeV/c 2 ) 



FIG. 3. The (a) KKirir and (b) i^7r mass spectrum from non-i^*X* events. In (b), the dashed 
line shows the k fitted to Kg (1430) and the full line is derived from K + K~tt + it~ phase space plus 
a fit to K* + Kir phase space. 
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FIG. 4. The contribution of (a) 0~, (b) 0+, (c) 1+, (d) 2+, (e) 2~ and (f) 4+ from the slice fit 
to J/V> -» 7i^*K* 
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FIG. 5. Comparison of the mass spectrum of K*K* with (a) all, (b) 0~, (c) 2 + and (d) 2" 
contributions from the best fit; the dashed curve shows the 0~ component fitted to the broad 0" 



resonance in Ref. [14]. 




FIG. 6. (a) Comparison of J/ip — ► jK*K* data and fit for (a) cos 9 K + + cos 6^+ (both entries 
summed) and (b) x> ( c ) an< i (d) corresponding angular distributions for CP (full histograms), + 
(dashed), 2+ (dotted) and 2~ (dash-dot). 
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TABLES 

TABLE I. Branching ratios (BR) integrated over all masses, but excluding the r\ c . 



Channel 


BR 


BR{J/^j -► jK+K-tt+it-) 


(2.1 ±0.1 ±0.6) x 10~ 3 


BR(J/i> -► jK*K*) 


(4.0 ± 0.3 ± 1.3) x 10~ 3 


BR(J/^j -► 7 _ ) x 5i?(0~ -► K*K*) 


(2.3 ± 0.2 ± 0.7) x 10~ 3 


BR(J/^j -► 7/2(1950)) x B#(/ 2 (1950) -► K*K*) 


(0.7 ±0.1 ±0.2) x 10~ 3 


BR{J/^ -► 7 2") x Bi?(2~ -► K*K*) 


(0.9 ±0.1 ±0.3) x 10- 3 



15 



